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Abstract 

c/3 ' We study theoretically the disorder-induced smearing of the density of states 

in a two-dimensional electron system taking into account a spin-orbit term 

3 ' in the Hamiltonian of a free electron. We show that the characteristic en- 

— ^ ■ ergy scale for the smearing increases with increasing the spin-orbit coupling. 

O ' We also demonstrate that in the limit of a strong spin-orbit coupling the 

diagrams with self-intersections give a parametrically small contribution to 

K^ ' the self-energy. As a result, the coherent potential approximation becomes 

OO \ 

1^ ' asymptotically exact in this limit. The tail of the density of states has the 

(N 

energy scale which is much smaller than the magnitude of the smearing. We 

^^ ' find the shape of the tail using the instanton approach. 
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=2 ! PACS numbers: 73.20.Dx, 73.20.Hb, 74.20-z 
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It is well known how the random potential smears the band-edge in a 2D system. In the 
case of a white noise potential with a correlator 

< \/(r)\/(r') >= 7(5(r - r') (1) 

the characteristic energy scale for the smearing is 

E2D = 7^ • (2) 

Deep in the tail {E < ,\E\ ':^ E2d) the density of states (DOS) falls off exponentially 

/ I Pi \ 
p{E) oc expy-i^—J , (3) 

where the numerical factor ^ is approximately ^ ~ 5.8cl'0 . The form of the tail (3) follows 
from the instanton approach developed in Refs. ^,|| (see also the books ||,H). The prefactor 
in Eq. (|^), including the numerical coefficient , was derived in Ref. |^. In the intermediate 
region, E ~ E2D, the exact form of the DOS is unknown. Within the coherent potential 
approximation it was studied in Ref. |I|. The autors of Ref. ^ have also performed the 
approximate matching of the coherent potential result and the tail (^. 

Spin-orbit (SO) interaction modifies the energy spectrum of 2D electrons. The origin of 
this modification is either the absence of inversion symmetry in the bulkQIa or the assymetry 
of the confinement potential. In the latter case the SO interaction can be taken into account 
by adding the following term to the Hamiltonian of a free electrons, 

Hso = a{& X k)n, (4) 

where the components of & are the Pauli matrices, n||z is the normal to the 2D plane, a is 
the SO coupling constant, and k stands for the electron wave vector. The energy spectrum 
of the Hamiltonian 



2 / |ifc2 a{h + zky)^ 



2m 
consists of two branches 



"so 



H = ^— fc2 ^ ^, 



(5) 



\a{kx-iky) ^k j 



E^(k) = a|k| , E2(k) = - — + a|k| 



2m 



The coressponding eigenstates have the form 



,(1,2) 
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^'^''\v)=e^^^x\i 



kr, (1,2) 
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where the spinors Xk ' ^ire defined as 
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Here 0k is the azimutal angle of the wave vector k. The lower branch, -Ei(k), has a minimum 

at 



k = kn 



am 



n 



2 ' 
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with a depth 



A 



ma 



2r 



(10) 



In the absence of a disorder the densities of states corresponding to each branch have the 
form 



p'^\E) 



rra 



1 + E/A + 1 
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pf\E) 



m 



1 + E/A - 1 



(11) 



'1 + E/A ■ - ' ' 27ih yi + E/A 

It is seen that p[ (E) is ID-like, in the sense, that it diverges as (— |-E| + A)~^/^. The energy 
spectrum (^ and the densities of states (|ll|) are shown in Fig. 1. 

The relation between the disorder and the SO coupling is measured by a dimensionless 
parameter 

E2D 7 



K 



2A a2 



(12) 



It is clear, that if k ^ 1, then the spin-orbit term has a negligible effect on the DOS. In 
other words, in the limit of weak SO coupling the smearing is still determined by the energy 



scale E2D- In the present paper we study the opposite hmit of a strong SO couphng (or 
weak disorder), k <^ 1. Remarkably, in this case the DOS can be found exactly. 

Let us first determine the characteristic energy scale, Eid, for disorder-induced broaden- 
ing. Using the golden rule, the relaxation time for an electron with energy close to E = —A 
can be written as 

--lpf{E). (13) 

'Te 

Then Eid can be found from the condition Eid ~ Ti/te^jj , yielding 

E.D = ^iiaf/' . (14) 

We see that for k ^ 1 the new energy scale is much bigger than E2D but much smaller than 
the depth of the minimum: 

^- = Si = «"=A . (15) 

This last condition allows a strong simplification in the calculation of the DOS. Indeed, Eq. 
(p!5|) suggests that the states in the region of smearing are composed of plane waves with 
magnitudes of wave vectors close to k^ , 



|k| - fco ~ ^j2mEiDlh ~ K,ko < k^ . (16) 

If we rewrite the energy spectrum Ei (k) as 

i^l(k) = -A + ^(|k|-A;o)^ (17) 

2m 

then Eq. (p!6D allows to consider the second term as a small correction. The crucial obser- 
vation, which allows the calculation of the DOS, 

in the closed form, is that under the condition k ^ 1 the contribution of the diagrams 
with self-intersections to the self-energy, Sk(-E'), is much smaller than the contribution 



of diagrams without self-intersections. In other words, in the strong SO couphng hmit the 
coherent potential approximation becomes asymptotically exact. To illustrate this statement, 
consider two second-order diagrams for the self-energy shown in Fig. 2. The contribution 
of the diagram (a) without self-intersection to ImTi can be written as 



The contribution of the diagram (b) with self-intersection is correspondingly 

*(i),.(i)Ar,.*a),.a) \i.A^) ,.(i)v,.*(i),.(i) 

k) 






(20) 



Noting that the scalar products (Xk Ak' ) ^^^ equal to 



*(1),X1)^ 



(xrxiV) = cos|r!l^)e-'^^. (21) 



the integration over the angles 0ki and 0k2 in Eq. ([I9| ) can be easily performed. The 



main contribution to the integrals over absolute values ki and k2 comes from the regions 
1^1 ~ ^0 1 ^ ^0 ) 1^2 — ^0 1 ^ ^0 • Then, using ([T7|), the energy denominators can be simplified 
to 

E - i5;i(ki,2) = i^; + A - |-(|ki,2| - Kf. (22) 

2m 
As a result we get the following estimate for ImYP^ 

In contrast to Eq. (|1^), in the second diagram the condition that the magnitudes of k, 
ki and k2 are close to /cq restricts the integration over angles. Indeed, consider the last 
energy denominator, ( _E + A — ^(|ki + k2 — k| — fco)^ ) , in Eq. (|20|). It is easy to see that 
|ki + k2 — k| can be close to k^ only in three domains: i) |ki — k| ^ /cq , m) |k2 — k| ^ /cq , 
III) |ki + k2| ^ fco . The size of these domains is determined from the condition 



|ki + k2-k| -A;o~ V^|^ + A|/r. (24) 

For the case t), for example, this condition confines the angle 0ki within the interval 



'm\E + A\/{nko). (25) 

Then the estimate for /mE^^^ yields 

/mS(^)^7^f „"^^° X '^. (26) 

' W\E + A\J ^ ' 

Thus, we get the following estimate for the ratio of diagrams (a) and (b) 



ImJ:W - yME + A\/{nko). (27) 

In the region of broadening, IE' + A| ~ Eid , this ratio is of the order of k^^^ <^ 1 . More 
accurate estimate (see below) gives K^^^\n{l/K). 

Once the diagrams with self-intersections can be neglected, the summation of the re- 
maining series is straightforward and yields the following equation for the self-energy 

It is easy to see that IiriL]^{E) does not depend on k. Although an explicit dependence on 
(/)k is present in the numerator of the integrand, it disappears after the angular integration. 
Substituting for -E(ki) the expansion ([T7|) and performing the integration, we obtain for 
JmS 

, ^ ErDJ2'l^ 



O • 



2^/3 M,,D 

where the energy e is defined as 

e = E + A- ReE , (30) 

and the dimensionless function f{x) satisfies the algebraic equation 

X + \ fixY + x^ 



The function /(x) is shown in Fig. 3. It turns to zero at a; = —2 ^/^. In the vicinity of 
X = —2^^/^ it exhibits a square-root behaviour 



^(^)^^--5^©"V(S. 



25/631/2 / I 

/(^)^^I7^\/^ + ^- (32) 

This behavior is usual for the coherent potential approximation. Using (0), the density of 
states ( ITSD can be expressed through the function f{x) as follows 

nj 2'Kh^\Kj \EiD 

Clearly, the vanishing of the DOS at e = —2^^^^Ei£, is the consequence of neglecting the 
diagrams with self-intersections. Taking these diagrams into account leads to the smearing 
of this singularity and formation of the tail of the DOS. The fact that intersecting diagrams 
are relatively small indicates that the characteristic energy for this smearing should be much 
smaller than -Eid- Indeed, below we demonstrate, using the instanton approach, that the 
DOS in the tail has the form 

^^^^"^"K-i^..ln(l/k|))- ^''^ 

It is seen from (^) that the rate of the decay of the DOS in the tail is E2D ln(A/_E2_D) ^ -Eid- 
Note that at |£:| ~ A, Eq. (|3^) matches the result (^ for the zero SO coupling. This 
conclusion could be anticipated since at energies \e\ ^ A the density of states does not 
depend on the SO coupling and Eq. (Q) applies. 

Within the instanton approach the density of states is given by 

p{E) <x exp^-^ J dhmr)\'y (35) 

where the function $(r) is the solution of the nonlinear equation 

H^{r) - |$(r)p$(r) = E$(r) . (36) 

When the energy E is close to —A, the two-component wave function $(r) is modulated in 
space with a period 2TT/kQ. Then it is convenient to perform the Fourier transformation of 
Eq. (|36D. Substituting 



<l>(r) = |d2rA(k)xL'^e-'^^ (37) 

we obtain 



Since A(k) depends only on the absolute value of k, the angular integration in (|38| ) can be 
easily performed. Using (^if) we obtain 



A{k)[E^{\i) - E) = 71^ fdrrfdkikJdk2k2fdkshA{ki)A{k2)A{k3) 

X [Joikr)Joik,r) + Mkr)J,{k,r)][Joik2r)Joik,r) + J,{kr) J,{k,r)] , (39) 

where Jo{x) and Ji(x) are the Bessel functions of the zeroth and the first order respectively. 
Now we make use of the fact that for E ^ —A the magnitudes of all wave vectors in 
Eq. (p9|) are close to fco- More precisely, the typical range of the change of each k is 



\k — ko\ r^ kir = Jm{E + A)/?i^. If we replace k, fci, k2, k^ in the arguments of Bessel 
functions by fco? then the product of Bessel functions will fall off as r~^ and the integral over 
r would diverge logarithmically at r — i> oo. This divergence should be cut at r ~ k~^. Then 
we obtain 

A{k){^^{k-k^f-e^ = AkoHko/k,)(^J^^dk'A{k')y . (40) 

The obvious solution of this equation is 

^{k-ko)^+\e\ 
Substituting (^) into (|40D, we find the value of the constant C 

^= 2^*""" (|?) ""'""'"<*»''*='■ (^2' 

Performing the inverse Fourier transformation, we obtain the solution of the instanton equa- 
tion in the coordinate space, which is valid for r ^ k^^: 
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X 1/2 f J-.(h.^\ \ 

$(r)=27r2Cfco' 



me\ 



Ji{kor) 
y -Jo(V) J 



(43) 



Finally, upon substituting (^) into ( pS]) we arrive at (Q). 

The above calculation allows to improve the estimate of the diagram (b) in Fig. 2. 
Note, that with A{k) in the form (^) the right-hand side of Eq. (^) coincides with the 
expression (^) for the diagram (b). This becomes obvious if the integration over r in 
Eq. (|38|) is performed explicitly. Comparing Eq. (|^) and Eq. (p6|), we see that in the 
order-of-magnitude estimate of the diagram (b) the factor ln(/co/^e) was missing. 

In conclusion, we have calculated the DOS for 2D electrons in the Gaussian random 
potential in the limit of a strong spin-orbit coupling. The summation of the diagram series 
became possible due to the fact that in the absence of disorder the energy spectrum has a 
minimum at some finite k = k^. This causes a singularity in the bare DOS. As a result, the 
magnitude of smearing of the DOS increases with increasing the SO coupling. Note, that 
the depth of the minimum in the energy spectrum decreases in the presence of a magnetic 
field and disappears completely when the Zeeman splitting exceeds 2 A. Then the smearing 
of the DOS occurs independently for both spin projections, so that the magnitude of the 
smearing becomes i?2D^the same as in the spinless case. 

We did not analyze the real part of the self-energy, i?eS. In fact, -ReS is determined 
by the large values of momentum and diverges logarithmically. This divergence is the con- 
sequence of the zero correlation radius of the random potential. For a small, but finite 
correlation radius -ReE causes a shift in the position of the band-edgeH'El. In the calculation 
of the shape of the tail this shift did not show up. This is a common situation for the 
instanton approach, within which the conversion from the bare to the "physical" energy 
occurs on the stage of calculating the prefactorH'ty. 

The applicability of the theory developed requires the SO-induced energy scale A to 
be bigger than the inverse relaxation time E2D in the absence of the SO coupling. This 



condition seems to be met in high mobility silicon MOSFETgi^tJ. According to Ref. O, the 



coupling constant a in this structures is ~ 2 ■ 10 ^niev-cm, which corresponds to A ~ IK. 
However, in the experimentally interesting situation where the metal-insulator transition 



occurs (see, e.g., recent references O-IT^), the Fermi energy lies much higher than A. 
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FIGURES 
FIG. 1. (a) The energy spectrum of a 2D system with spin-orbit couphng. (b) The density of 

states for two branches of the spectrum. 

FIG. 2. Two second-order diagrams for the self energy S. (a) The diagram without 
self-intersection, (b) The diagram with self-intersection. 

FIG. 3. Dimensionless function f{x) defined by Eq. (31). 



13 



a) 



^ 



N 



W 



-A 




Fig. 1 



b) 



o 

Q. 




Fig. 1 



a) 



/ 

/ /' 

/ / 

I I 



k k. 



kj k 



b) 



\ \ 

\ \ 

\ \ 

I I 



k kj kj + kg-k kg k 



Fig. 2 



1.25 



0.75 



X 



0.5 



0.25 








-2 



4 
X 



8 10 



Fig. 3 



